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between metastable and all possible configurations, for both ferromagnetic and bond-disordered 
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I. INTRODUCTION 

Small-world systems exhibit remarkable cooperation, 
not found in complex systems with e.g. an Erdos-Renyi 
structure. The origins of the appellation 'small-world' 
can be traced to the now famous experiment by the 
Harvard social psychologist Stanley Milgram 0. The 
outcome of this experiment pointed to the fact that the 
structure of many real networks is such that distant nodes 
can in fact be connected via long-range short-cuts. This 
architectural property leads to small path lengths be- 
tween any pair of nodes and thus enhances information 
processing and cooperation. 

The importance and ubiquitous nature of small-world 
structures in complex networks received further attention 
by the seminal paper of Watts and Strogatz 0] in which 
the authors proposed the small-world structure as a way 
to interpolate between so-called 'regular' and 'random' 
networks. Surprisingly, the simple small- world architec- 
ture can be found in many different circumstances, rang- 
ing from linguistic, epidemic and social networks to the 
world- wide- web. By now, a large body of work has been 
devoted to the study of small- world networks, mainly nu- 
merical, with emphasis e.g. on biophysical 0,0,0 or so- 
cial networks 0] and to a lesser extent analytical • 
For recent reviews in the area of small-worlds see e.g. the 
articles 0113, or the books [HG2. 

From a statistical mechanical point of view such sys- 
tems combine two universality classes: a sparse 'graph' 
structure, which is superimposed upon a one-dimensional 
'ring'. Thus, every node on the ring has a local neighbor- 
hood and a certain number of long-range connections to 
distant parts of the chain. It was shown in (and also in 
[l3l | for the case of XY spins) that this construction sig- 
nificantly enlarges the region in parameter space where 
ferromagnetism occurs. In particular, it was shown that 
the ferromagnetic-paramagnetic transition always occurs 
at a finite temperature for any value of the average long- 
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range connectivity (however small). On a technical level, 
in evaluating the relevant disorder-averaged free energy 
one is immediately confronted with the problem of diago- 
nalising a 2" x 2 n transfer matrix where n represents the 
replica dimension. Although for an infinite system size 
obtaining only the largest eigenvalue suffices for the eval- 
uation of the free ene rgy, one can in principle follow the 
systematic analysis of [l4j to derive the entire spectrum of 
eigenvalues and thus evaluate e.g. correlation functions. 

Several important issues remain to be understood for 
small- world systems. In this paper we evaluate the num- 
ber of metastable configurations, or more precisely, the 
number of equilibrium configurations in which spins align 
to their local fields. With this definition the energy of 
the system in a metastable state cannot be decreased by 
flipping a single spin. Such configurations can be e.g. re- 
sponsible for trapping the microscopic update dynamics 
in locally stable states. Thus, from an experimental point 
of view it is advantageous to know what is the relevant 
size in phase space occupied by such states. 

The computation of the number of metastable con- 
figurations in small-world systems is generally an in- 
volved problem, both analytically and numerically In- 
dicatively, on sparse random graphs structures without 
the superimposed 'ring' the evaluation of the number 
of metastable configurations, or the so-called configura- 
tional entropy, has only recently taken off 0, Il6l 
following the course of the relevant analytic techniques 
(as e.g. in 0,0, H(j)- In particular, using the replica 
method, the solution of the ferromagnetic Poisson graph 
has been studied in while with the cavity method 
the authors of ^(| examined the bond-disordered Bcthe 
lattice. These results agree well with results of numeri- 
cal enumerations 0, |22j and also serve as good limiting 
tests of our findings. 

This paper is organized as follows: in Scction[n]we first 
define the small- world model. In Section IlIII we express 
the generating function of the system as a saddle-point 
problem which we then evaluate in Section llVl using repli- 
cated transfer matrix techniques and within the replica 
symmetric approximation. In Section [V] we present our 
results and finally we present some concluding remarks 
in Section IVTl 
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II. THE MODEL 

Our model describes a system of N Ising spins er = 
(a"x, . . . , ctn)-> with cTj G { — 1,1}, arranged on a one- 
dimensional lattice. There are two different couplings in 
this system: firstly, nearest-neighbor interactions of uni- 
form strength J and secondly, sparse long-range ones. 
To model the latter we will assign the random variable 
Cij for every pair of sites representing whether a con- 
nection exists (cy = 1) or not (cy = 0). This variable 
will be taken for all i < j from the distribution 
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so that, on average, every site has c long-range connec- 
tions. In the small-world context one takes c to be a small 
number of order 0(1) while c/N — ► 0. The bond-strength 
Jij of the long-range coupling between any pair of spins 
(with i < j) will be taken from the distribution 

QAJii)=pSj tj ,j + (l-p)6 J(j ,-j, (2) 

for some J > 0, so that p = 1 corresponds to a model with 
strictly ferromagnetic interactions. To allow for detailed 
balance we will also consider absence of self-interactions 
and symmetry of the connectivity matrix, namely cu = 0, 
Cij = Cji and = Jj\. At thermal equilibrium the above 
system can be described by the Hamiltonian: 



ft(<r) = -i^cr 4 /i 4 (er), 



(3) 



with the local fields defined as: 
hi(cr) = ^2 \ Jo(8j,i+i + 



-J,, 



We now impose the condition for metastability: in a sim- 
ilar spirit as e.g. in [23l l24j we call a configuration er 
mctastablc if all spins align to their local fields, i.e. 
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where 0(x) = 1 for x > and O(x) = otherwise. 
Notice that we have taken O(0) = 1 which is dictated 
by the underlying physics: for spins which receive a zero 
local field, the energetic cost of aligning to either of the 
two possible directions is identical. A consequence of 
the above definition is that any metastable configuration 
er MS is a local minimum of the Hamiltonian: 
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Since we are interested in evaluating the number of 
mctastablc configurations we will define the following 
generating function: 



N 
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where (3 represents the inverse temperature and c = 
{0^}, J = {Jij}- As in [I5I Ht|. the entropy density s m 
of the metastable configurations can be evaluated from 
61) via a m =pdp{pf)-pf. 



III. SADDLE-POINT EQUATIONS 

To evaluate the disorder average in (0J we begin by in- 
voking the replica identity: (logZ) = lim n ^o ^l°g(Z n ). 
As the disorder variables {cy , Jij } lie within the 
©-function we insert the following unity into our expres- 
sion 



1 = jl[dhf 5[K-hi(<r a )}, 



(5) 



with er" = (erf, • • • , afj) and a = 1, . . . , n. This allows 
us to conveniently re-locate cy , Jy into exponents where 
averages can be taken more easily, namely 
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where {dhf dhf } = J], a {2ir)- 1 dhf dhf . Let us concen- 
trate on the last line of the above which contains the 
disorder. After symmetrizing with respect to the sites 
i < j it leads for A — > 00 to: 
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where (■ ■ ■} j denotes an average over the binary random 
variable J taken from the distribution Qj(-), eq. J5J). 
We have used the fact that c/N — ► to recast the result 
of averaging over {c^} into an exponential form. Note 
now that upon inserting the unities 1 = Scr cri and 
1 = ^2i T St.cTj where o, t are auxiliary vectors in replica 
space (and we have denoted er^ = (<jj, . . . , erf)), one has 
effectively created an order function Pj((t,t). As usual 
it can be inserted into our generating function via: 



i=/nn dPj(cr,T) 
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As in [23, to understand the physical meaning of the 
above order function one needs to add a generating term 
in the replicated Hamiltonian ^ Q H(cr a ) — > J2 a H(er a )+ 
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riPj(cr, t) and take the derivative df /drj\ v= o in Q. One 
then sees that upon introducing the identities © the or- 
der function becomes the distribution of replicated spins 
with one connection removed (equivalcntly, it becomes 
the distribution of replicated 'cavity' spins). 

We now aim to eliminate from our expressions the set 



of fields {hf, hf }. This can be done by replacing the delta 
function in JHJ by its Fourier representation (for details 
see appendix 0). As an end result we obtain an extrem- 
isation problem over the density Pj(ct,t) expressed in 
terms of a trace over a transfer function, namely 
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where J G { — J, J} and with the abbreviation 
tr(T N [P]) = tr (T n [Pj,P_j}). The order functions P ±J 
are to be evaluated self-consistcntly from 



j = l (J i— <T N 

x F^_ ua . t(T . +1 [P, a, t] [] /rr a ,7 [P] (10) 

We have absorbed the normalisation constant 
lli r o' i _ 11 <7- i ,(T i+1 [-P] in the proportionality 
symbol. The traces in (JSJ and (|10|l involve correlations 
between next-nearest neighbors and can be evaluated in 
a spirit similar to the transfer matrix technique. The 
relevant tensor is defined over a 2 n x 2" x 2™ space with 
elements 
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with the convention that the fi = term of the above 
equals e~ c ]^[ a Rf(0). The quantity F that specifies our 
order function is 
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We have used the abbreviation 
Rf(x) = 



(12) 
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e |/3^(^oK ? +1 +^- 1 )+S) [ (T «(J o ( (7 « +1 + + -)]. 

Note that due to our symmctrization with respect to 
site indices in (JJJ we have ended up with a symmetric 
quantity, namely Ti, Xi = T Xi where L is the 3x3 matrix 
Lik = S e+ k,4 and Xi = (<x,_i, <Xj, cr i+1 ). 



The structure of the function Rf (x) indicates that the 
input x is related to the long-range field received by a 
site i. The tensors T and F differ only in their input 
to this function. Since this is proportional to ±J/c we 
understand that it is related to the 'effective' (or, 'cavity') 
and true local-field, respectively. 



IV. REPLICATED TRANSFER MATRIX 
ANALYSIS 

To interpret the spin summations as matrix multiplica- 
tions in equations @ and IjlUfl we need to transform our 
variables such that the traces in these equations involve 
nearest-neighbor correlations only. This can be done in 
several ways. For instance, let us introduce the auxiliary 
spins: 



/ a.l a,2\ 
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To suppress the replica index above, we will occasion- 
ally use the more compact notation s.; = (s^ 1 , s\ ) = 
(<Ti, <Ji+i) with vectors now defined in replica space. 
With the above, we can now transform T (Ti _ 1 .cr i ,(T 1+1 into 
Ts i _ 1 ^a i and in particular 
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with e = (1, 0). Similarly for F: 
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The Kronecker deltas on the right-hand side of the above 
expressions impose the transformation l|14|) . Similarly to 
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equation l|13|l we now have 



af ■ (Josti + -e) 



This transfer matrix is also symmetric. In particu- 
lar, it obeys T^y = Ty where L is the 4x4 ma- 
trix Lik — (5^+^.5 (the so-called Dirac matrix En) and 
y i = (sj,Sj_|_i). With the above definitions we can now 
write the self-consistent equation (|lt)(l in a transparent 
way: 



ti (Q^[P, a', tT,r] T"- 1 ^] 



tr (T N [P}) 



(16) 



with the auxiliary matrix 
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(17) 

To proceed with the evaluation of the traces involved in 
(|9"|) and l|16(l we now aim at diagonalizing the transfer 
matrix T. Our analysis hereafter will follow closely @- 
To this end, let us consider the eigenvector equation cor- 
responding to the largest eigenvalue Ao(n), namely 



J2Tss-[P] u(s') 

S' 



A (n) u(s). 



(18) 



Note that we have only defined a 'right' eigenvector. It is 
sufficient due the symmetry of our transfer matrix. Next, 
using l|17|l and i|18|) wc can rewrite the self-consistent 
equation i|16|) in terms of it(s) 



(19) 
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Equations (|18|) and l|19|) are the basis of our analysis in 
the subsequent sections. 



A. Replica symmetry and self-consistent equations 

Since the order function P depends on n via the dimen- 
sionality of its arguments we must now make an ansatz 
that will allow us eventually to take the limit n — ► 0. The 
simplest choice corresponds to considering permutation 
invariance of P with respect to its replica indices. This 
symmetry is guaranteed by considering e.g. the following 
form: 

dw ^ h) wm> ' (20) 

with the short-hand notation h = (hi, h,2, hs), dWj(h) = 
dhWj{h) and Af(h) ensures that J dhWj(h) = 1. We 
also assume that the eigenvector u(s) takes the form 



u(s) = / d$(x) e 
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(21) 



With these assumptions we can now proceed further and 
rewrite the extrcmisation problem in terms of the pair 
of densities $ and Wj. The starting point is equations 
(|18fl and (|19|) respectively. Inserting our assumptions 
(|20|l and (|21|l leads after some algebra to the following 
set of closed equations for n — > 0: 
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In the above we have used the function fi (a, r) with 

A (c, r) = a, f 2 (a, t) = t, f 3 (a, r) = err, 
while G, H correspond to 
G„, T (x,h) = 

lo g| E E e^+^+^Sr^iJ^ J»r v ) 1(24) 
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Tl — Tp r/,ui=± 

V<ft I 



and 
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Finally, to calculate © we need to determine the 
largest eigenvalue \o(n) for n — > 0. The starting point 
here is our eigenvector equation (|18J) . Evaluating the 
traces over the spin variables with the definitions of the 
transfer matrix l|15|) and eigenvectors (|21|) leads for n — > 
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to 
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so that Ao(0) = 1. We have defined 



a.cr 1 ' — ± 

tf(a;, = J] H e /3(xi<Tl+X 2 a3+ S3 <71<73) 

^CTl ,CJ2 ,C3 ( 1 ^ ) ■ 

With the expression (|26[) we can now evaluate J!5J|: 



(28) 



dWXh) dWj{ti) logD(h, ti 



J 



/ e c ' 



J da;' $(a?') J 



K(x,{K}) 



A(x) 



(29) 
■••Sf* 



with 
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a,a' — ± 

which is our final result. From this equation we can in- 
spect the physical meaning of the densities Wj (h) and 
<$>(x). The Poisson distribution e _c c M //i! of mean c can 
be clearly associated to the degree distribution of the 
graph. Once a degree has been sampled from this distri- 
bution one performs /1 integrals over the densities {W jj\ 
and one over <E>(a;). Thus, we can think of the Wj v (h) 
as the distribution of 'effective' fields (or, so-called 'mes- 
sages') coming from the long-range connections and $(a;) 
as those coming from the ring neighborhood. 



B. Benchmark tests of the theory 

Given the complicated structure of our equations, we 
now wish to inspect the validity of the theory against 
simple benchmark tests. 

Firstly, in the absence of the 'ring' structure, and for 
strictly ferromagnetic interactions, the equations must 
reduce to those found in Ref. ^5|. Indeed, setting into 
the update functions I|24I25|I Jo = and p = 1 we find 
that after the first iteration of (|22f) the density $ collapses 



to X2, X3) = 8(xi)4>(x2)8{xz). Details of the non- 

trivial function <fi are not important for the purposes of 
this section. Filling in this information in the right-hand 
side of l)23|l leads to several simplifications as a result of 
which the dependence of Wj(h) on (j>(x) drops out com- 
pletely. The resulting closed equation is the one found 
in 15]. Thus, at the level of the self-consistent equation 
the expressions reproduce the correct result. Next, we 
consider the free energy. Clearly, the energetic part of 
(|%9l depends only on Wj{h) and in the special bench- 
mark case takes the same form as the energetic term of 
|l5l| . The entropic term on the other hand, depends ex- 
plicitly on the reduced density of fields 4>(x) which is 
coupled to the functions A(x) and K(x,h), eq. I|27|) and 
(12811 respectively. Here, it turns out that one can write 
K(x, h) = A(x) K{h) which effectively removes A{x) and 
(f>(x) completely from (I29|l . The resulting expression re- 
produces the free energy of [l5j . 

A second test of the theory is again st the small-world 
thermodynamic analysis of Ref. |8|, l25l | . To map the gen- 
erating function |0J to the free energy of that system 
we set O(x) = 1 for all x. This removes the stabil- 
ity condition from our definitions. After the first iter- 
ation of (|22fl we now find that the function $(a;) col- 
lapses to $(xi, X2, X3) = 8(xi)<fi(x2)S(x3 — \ Jo) and us- 
ing this to iterate we obtain that W(hi,h 2 ,h 3 ) = 
w(hi)5(h 2 )S(h 3 — J/2c). Thus, in both cases only one of 
the three components is non-trivially distributed. With 
these relations we recover at the second iteration of 
(|22I23[) the self-consistent equations of [25|. Equations 
(I22I23[> also reduce to (29,30,32) of if the analysis of 
[K would have been based on symmetric transfer ma- 
trices. In this case the final result l|29[) reduces to the 
correct free energy. 

Finally, by inspection of the physical interpretation of 
equations (1221231) and (12911 we can map our mode 1 ! onto 
the one of [T(| which evaluates the number of metastable 
configurations on a Bethe lattice. This can be done by 
appropriately converting the Poisson degree distribution 
to a Bcthe-lattice one. We have done this test numer- 
ically and within the limits of precision we find good 
agreement with the results of fig . 



V. RESULTS 

We are now interested in obtaining the energy e and 
entropy densities s of the metastable states. These 
can be generated from / through simple relations, i.e. 
e = dp{(5f) and s = (i(e — fY. To obtain these we have 
solved equations (|22|l and i|23|) through 'population dy- 
namics' [l8| and used simple Monte Carlo integration 
recipes to evaluate (|29fl (typically, in executing popula- 
tion dynamics the population size has been taken of the 
order of 10 5 and we assumed algorithmic equilibration 
after 100 steps). Since the profiles of Wj(h) and <fr(x) 
depend on the temperature, differentiation of / with re- 
spect to (3 will involve derivatives of both of these den- 



sities. One of these, namely dwf • dpW, trivially van- 
ishes as we are at saddle-points of the order function 
Pj(ct,t) and, consequently, also of Wj(h). The deriva- 
tive d$f ■ 9/3$ however, may not necessarily vanish as 
we have not extremized / with respect to Therefore, 
we cannot assume that the energy e = 9^(/3/) is given 
by a simple partial differentiation of l|29|) (which indeed 
leads to incorrect results). To proceed analytically, one 
is required to obtain further closed relations for <3> = d/3<$> 
which, given the complexity of the equations involved, 
appears to be a hard task. Here, we have chosen to carry 
out the differentiations numerically. In all cases we have 
taken the average connectivity to be c = 2. 

Let us now describe the results. First we take the sim- 
plest case where long-range interactions are of uniform 
strength, i.e. p = 1. Thus, the only source of disorder 
in the system comes through the connectivity variables 
{c.jj}. In figure ^ we plot the energy against inverse tem- 
perature for three different values of Jo and with J = 1. 
In each case we compare the energy density of the system 
e when we allow all configurations to be visited in phase 
space (i.e. with 0(x) = 1 for all x) against the energy 
e m of only the metastable subset of configurations (with 
Q(x) = 1 if x > and otherwise). 

Since, by definition, metastable states are minima of 
the energy landscape we expect that e m < e at any given 
temperature which indeed is verified by the numerics. As 
we increase the strength of the short-range couplings Jo 
the system will typically require a higher noise level to 
destroy the order. For this reason we see that the lo- 
cation of the phase transition towards low energy values 
decreases with (3 as Jo is increased. For /3 — * oo one can 
find the ground state energy of the system by simple in- 
spection of the Hamiltonian (0, namely e gr = — Jo — J/ 2. 
Furthcmore, since in the regime of low temperatures we 
expect the system to be in a locally stable state we can 
anticipate that e ~ e m . These physical arguments are 
in agreement with the numerical results of figure Q On 
the other hand, for (3 — > 0, noise dominates the micro- 
scopic spin dynamics and thus the energy of the system 
e typically averages to zero for all values of Jo- We also 
observe that the transition to the ordered phase is less 
smooth for e ?n than for e. This effect which has also been 
reported in for the special case of Jo = is due to 
non-linearities induced by the Heaviside function. 

In figure |21 we plot the entropy against the energy for 
different values of Jo, and with p = 1, J = 1, c = 2. The 
low-energy part of the graph corresponds to regimes of 
low temperatures. As before, we compare the entropy s 
that would follow from a thermodynamic calculation in 
the entire configuration space against the entropy s m of 
the metastable states. First, we see that one always has 
s m < s, as one would expect. The difference between 
the two entropies varies significantly across the energy 
axis. For instance, for high energy values (where the tem- 
perature is practically infinite) this difference reaches its 
maximum value. On the other hand, for low tempera- 
tures, both entropies reach their minimum value which 





\ ■ . 

\ 

\ 

— \ 

"V-.. 






\ 

\ 

\ 

\ 


— Jo = 
■ ■ ■ Jo = 1 

- - Jo = 2 










\ \ 

\ \ 
e 1 \ 

Cm , ^ 

1,1,1 







p 

FIG. 1: Comparison of the energy e of a small- world sys- 
tem in which spin configurations can visit all possible config- 
urations against the energy e m of the metastable configura- 
tions, plotted as function of the inverse temperature (3 and 
for Jo = {0, 1, 2}. We see that e m < e and that for f) — » oo 
e m ~ e. Parameter values: J = 1, p = 1 and c = 2. 
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FIG. 2: Comparison of the entropy s of a small- world system 
in which spin configurations can visit all possible configura- 
tions against the entropy s m of the metastable configurations, 
plotted against their energy e,e m , for Jo = {0,1,2}. For 
Jo = the ground state entropy has a finite value which van- 
ishes as soon as Jo > 0. Parameters values: J = 1, p = 1 and 
c = 2. 



for any Jo > is zero. In the special case where Jo = 
the graph will typically consist of disconnected clusters 
which causes the observed degeneracy. However, as soon 
as the 'ring' connects all spins together this degeneracy 
is lost and the ground state entropy is zero. 

Let us now examine the case of bond-disorder. In figure 
|3]we present the entropies s and s m against the energy 
for different values of Jo and with J = l, c = 2, p = 
1/2. Firstly, we observe that the ground state energy is 
significantly higher compared to the case of p = 1 . This is 
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FIG. 3: We show the entropies s and s m for small- worlds with 
bond-disorder. In this case the degeneracy of the ground state 
remains finite even for Jo > 0. Parameter values: J = 1, 
p = 1/2 and c = 2. 
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and c = 2. 



fraction of sites with a zero local field. However, as one 
increases the strength of the (ferromagnetic) short-range 
couplings this fraction of sites becomes smaller and the 
degeneracy of the ground state gradually disappears. To 
illustrate this effect we plot in figure 0] the ground state 
entropy against the short-range coupling strength Jo. We 
also see a 'jump' precisely at Jo = 0. 

VI. DISCUSSION 

In recent years, the theory of complex networks has 
witnessed a remarkable growth. Within the area of com- 
plex systems, the special subset of 'small-worlds', has 
aroused the curiosity of theorists and experimentalists 
alike due to the striking coopcrativity phenomena that 
they allow. In particular, for any value of the average 
long-range connectivity (however small), small- world net- 
works can have a phase transition to an ordered phase at 
a finite temperature. Small world architectures have been 
observed in a wide range of real complex systems. 

For a theorist, several important questions arise re- 
garding the emergent collective properties on such sys- 
tems. In this paper we have evaluated the number 
of metastable configurations. In a spirit similar to 
D3 IHl the 'metastability' condition constrains the 
partition sum over configurations in which spins align to 
their local fields. From an analytic point of view, two 
are the main stumbling blocks: firstly, the non- linear na- 
ture of the stability condition, and secondly, the diago- 
nalisation of the relevant transfer matrix. Our numeri- 
cal results suggest that, for low temperatures and in the 
case of bond-disorder, the metastable configurations tend 
to dominate the space of equilibrium states. We also 
see that superimposing the one-dimensional 'backbone' 
structure leads to a significantly smaller degeneracy of 
the ground state (which, in fact, vanishes for strictly fer- 
romagnetic couplings). As function of the short-range 
coupling Jo, there is a jump in the ground state entropy 
exactly at Jo = which is due to the formation of dis- 
connected clusters within the graph. 



due to the value of the local fields which will on average 
be smaller for p < 1 than for p = 1. We also observe 
that the ground state entropy can take a non-zero finite 
value even at Jo > 0. This is due to the presence of anti- 
ferromagnetic couplings in the system which increases the 
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APPENDIX A: DERIVATION OF THE 
SADDLE-POINT EXPRESSION © 

Our starting point are equations l|bl8[) . We replace the 
delta function in JSJ) by its integral representation which 
results in 



c,J 



A II dPj(<T,T)dPj(<T,T)] C Xp 
<T,T,J 

i —~ J a 
{dhf} 



l er.T j Vctt / 



N 

n 

i=i 



(27t)™ 



li>0 ^ Ti—TpJi—Jp V</i 



We now see that the last line of the above expression variables {hf} can be done immediately. The result is a 
has factorised over site indices and the integral over the delta function which we use to eliminate {hf}: 



9 



( Zn )c,J = fill dPj(*,T)dPj(cr,T)}exp *^^P t7 ( ( r,r)F i7 (,T,r) + ^(^(P J ( ( T,r)P J (r,,T)) :r ) 

J (T,T,J L (T,T J \(TT / 



N 



x / I] e^^^eCof ft?) 5 ft? - J (a? +1 + o?^) - i £ Xr« 1 

a L v=\ J J 

= /[II ^(^T^^.Tjlexp i^^P >7 ( < T ) T)P i7 ( e 7,T) + ^(^(P >7 ( -,T)P >7 (T, ( T)) 

cr,T,j L ct,t j \<tt 



X 



E IHE^r E E n^>— )] 



With the above, the generating function Q can be writ- 
ten as an extremisation problem, namely: 



-Pf = lim -Extr 



£ £ P^, r)P J ( < r, r) + |X! < P ^ CT ' T ^ T > + ^ ^ l0 § tr ( f ^) 



0",T J 



with the transfer function T[P] given by 



(Al) 



(A2) 



T, 



p \ = E^ E E m^(^)] 



]J e §/3<(M<^+<-iM£,=i^O e (af(Jo(a? +1 +a?_ x ) + i J^r?) j ( A3 ) 



Variation of l|A2() with respect to the function Pj gives {—J, J}. Using this identity in (|A2|) leads to (©• 
the relation Pj{a 1 r) — icQj(J) Pj(t,ct) for J = 



